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In this paper, a new methodology for Transmission Expansion Planning (TEP) in deregulated electricity 
market is presented. The proposed TEP is associated with Reactive Power Planning (RPP), reliability 
assessment and also consideration of wind and load uncertainties. The proposed planning aims at 
investment cost minimization, social welfare maximization and satisfying reliability constraint at the 
same time with taking into account wind and load uncertainties. Expected Energy Not Supplied (EENS) is 
used as an index for reliability evaluation. At first, Monte-Carlo simulation is used to obtain the 
Probability Density Function (PDF) of Wind Turbine Generator (WTG) output. Then, the WTG and load 
uncertainties are considered in TEP formulation. Particle Swarm Optimization (PSO) method is 
considered to solve the proposed planning problem which is a constrained nonlinear mixed integer 
optimization programming. Simulation results on two standard test systems (Garver and RTS systems) 
verify the effectiveness of the proposed planning for consideration of wind and load uncertainties in TEP 
problem under electricity market environment. Also, the proposed method leads to reduction of the total 
investment cost, the reliability improvement and the social welfare maximization. 

© 2013 Elsevier Ltd. All rights reserved. 
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Nomenclature 


Symbols and Indexes 


bi shunt susceptance of line or transformer ij (if ij is a 
transformer bj =0 )(p.u.) 

p” shunt susceptance at bus i (p.u.) 

Cok and Cık installation costs and unit costs for a VAr-plant at 
bus k ($) 


cı and c learning factors 

fi(q.u) cost function of locally reactive sources ($) 

fo(Ps, PI, IIs) cost function of social benefit ($/h) 

LZoest(t) global best value 

gj and byconductance and susceptance of the transmission line 
or transformer ij (p.u.) 


i,j bus indices 

Ng set of all buses 

nj number of new added transmission lines to branch j 
me maximum number of new added lines to each branch 
N, number of scenarios 

Ne project life-time 

Phes(t) i™ particle best solution 


position of the d‘ dimension of the ith particle in tth 

iteration 

pe® and Q% maximum limit of real and reactive power 
generation limits (p.u.) 

PY and Q™* minimum limit of real and reactive power 

generation limits (p.u.) 


Pia(t) 


q total amount of locally reactive sources (MVAr) 
qk amount of locally reactive source in bus k (MVAr) 
q™™ and q™™" maximum and minimum amounts of reactive 


sources (p.u.) 

Tt discount rate 

sfem and S° apparent power flow through the branches in both 
terminals (p.u.) 

Se apparent power flow limits (MVA) 


Vo total investment due to the addition of new circuits ($) 

vı total investment of locally reactive sources ($) 

v2 total cost of social welfare ($/h) 

Vannual annual cost of the proposed planning problem ($/year) 

Vid(t) velocity of the dth dimension of the ith particle in tth 
iteration 

ymax and v™" maximum and minimum voltage magnitudes 
(p.u.) 

Vtotal total cost of the proposed planning problem ($) 

w inertia factor 


Q set of all load buses 
keQ, kth load bus 


Vectors and matrixes 


c lines cost vector ($) 

n added lines vector 

N matrices containing the new lines 

No matrices containing the existing lines 

Pp and Qp real and reactive power demand vectors (MW 
and MVAr) 

Pc and Qg real and reactive power generation vectors (MW 
and MVAr) 

Ps; and Py; vectors of supply and demand powers in bus 
i (MWh) 

u binary vector that indicates whether or not to install 
reactive power sources 

Uk binary vector that indicates whether or not to install 


reactive power sources at bus k 

Vand © magnitude and angle of voltages vectors (p.u. and 
radian) 

Is; and I; vectors of supply and demand bids in bus i 
($/MWh) 


Abbreviations 


DG distributed generation 

EENS expected energy not supplied 

EENS™** maximum amount of expected energy not supplied 
FOR forced outage ratio 


GA genetic algorithms 

ISO independent system operator 
LMP local marginal prices 

LOLE loss of load expectation 


LOLE™™ maximum amount of loss of load expectation 


MCS Monte-Carlo simulation 

MTTR minimum time to repair 

OPF optimal power flow 

PDF probability distribution function 
PSO particle swarm optimization 

RPP reactive power planning 

RTS reliability test system 

TEP transmission expansion planning 
WTG wind turbine generator 


1. Introduction 


TEP problem aims at expanding the power system transmission 
network to serve the growing demand in the future. The TEP problem 
denotes where and when, new lines should be installed in the 
network to support the customers demand. The TEP problem in 
regulated power systems intends to minimize the investment cost, but 
in the deregulated electricity markets, the TEP plans to provide a 
competitive environment for all participants. 

TEP problem is a nonlinear mixed integer constrained pro- 
gramming. The commonly used models for TEP problem modeling 
are DC [1] and AC models [2,3]. Also, many different methods have 
been applied to solve the TEP problem. These approaches are 
divided into mathematical and Meta-heuristic approaches. The 
mathematical methods such as Linear programming [4], nonlinear 
programming [5], mixed integer programming [6] and Bender's 
decomposition |7] have been carried out for solving the TEP 


problem. Also, some heuristic optimization methods such as PSO 
[8], Tabu Search [9], Harmony Search | 10], Genetic Algorithms [11], 
Decimal Coded Genetic Algorithms [12], Genetic Algorithms-based 
quadratic programming [13], Chaotic [14], Ant Colony [15] and 
Differential Evolution [16] have been successfully used for solving 
the TEP problem. 

From another view, in TEP problem, many different parameters 
have been considered as objective function. The minimization of 
the investment cost is a conventional objective function for TEP 
problem [2]. Transmission surplus capacity as an objective func- 
tion for TEP problem has been investigated in [14]. The maximiza- 
tion of the transmission reliability and the minimization of the 
investment cost are considered as objective functions for TEP 
problem in [15]. Combination of TEP problem with RPP is pre- 
sented in [17]. This paper has shown the TEP problem associated 
with RPP results in fewer new transmission lines to be installed in 
comparison with the conventional TEP method. 
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In deregulated electricity market the objective functions are 
different from the conventional markets. In deregulated markets, 
the market concepts such as congestion surplus, social welfare and 
nodal prices are considered in the TEP problem. Considering the 
investment cost and the annual generation cost in the TEP 
problem under electricity market has been reported in [16]. In 
[17] the combination of TEP and RPP has been presented in 
electricity market. Coordination of transmission and generation 
capacity planning in a deregulated market has been investigated in 
[18]. A bi-level method to TEP under market environment has been 
studied in [19]. A market oriented TEP problem with considering 
the social welfare and investment costs has been reported in [4]. 

Recently, DGs such as wind farms are rapidly developed in 
power systems, because of their favorable characteristics [20]. 
However, wind power output is continuously changing and thus, 
an extra factor of uncertainties is introduced in power system 
operation and planning. Also, the loads in power system are 
uncertain. Therefore, in power systems with wind power and load 
uncertainties, the deterministic TEP methods are not suitable. A 
method with considering the wind power and load uncertainties 
should be incorporated [21]. 

In this paper, the AC-TEP associated with RPP is considered. 
The proposed planning is carried out in a deregulated electricity 
market with taking into account wind and load uncertainties. The 
proposed planning aims at the investment cost minimization and 
the social welfare maximization by considering reliability index at 
the same time. The reliability of planning is evaluated by using 
EENS index, which is calculated by using Monte-Carlo simulation. 
Also, the wind and load uncertainties are incorporated by using 
Monte-Carlo simulation. Also, PSO algorithm is used to solve the 
problem. Simulation results, which are carried out on two typical 
power systems, show the validity of the proposed planning. 

It should be noted that load uncertainties are modeled as 
normal Probability Distribution Function (PDF). The normal PDF 
is the most commonly used method to model different uncertain- 
ties in electric power systems such as load, cost, fuel etc.; it is also 
used in this paper [17-18,22-24]. 


2. Probabilistic WTG output model 


Wind turbine output is nonlinearly related to the wind speed. 
With changing wind speed, WTG output may vary between zero to 
its rated output and, hence leads to fluctuations in power flow. As 
a result, all power system characteristics such as nodal prices, 
social welfare, line transmission powers change. Therefore, these 
uncertainties should be incorporated in power flow formulation. 
A probabilistic model is a very suitable model for tackling this 
problem. The commonly used model for wind speed is Weibull 
distributions, and the shape and scale parameters of the distribu- 
tions can be derived from the mean and standard deviation of the 
wind speed [21]. The nonlinear model of a WTG relating the 
power output characteristics to the input wind speed has three 


0 Vi. V 


cl rate 


Fig. 1. Relationship between wind speed and WTG output power. 


parameters namely cut in speed (V,;), rated speed (Vate) and cut 
out speed (Vo). The output power of a WTG unit may be 
approximated by (1). The relationship between wind speed and 
WTG active power output is depicted in Fig. 1 [21]. 


(0) O<V<Vaä 

P= Prate P> Va <V< Vrate (1) 
Prate Vrate < V < Veo 
0 Veo <V 


By using wind speed distribution and WTG characteristic, MCS can 
be applied to simulate the WTG output distribution. For example, 
for a 150 MW WTG unit, the WTG output probabilistic distribution 
is depicted in Fig. 2. In this figure, the WTG parameters are as 
Vai=4, Vrate=10 and V,.=22 m/s, wind speed mean 5.4 m/s 
and wind speed standard deviation 2.7 m/s [21]. The WTG 
output distribution in Fig. 2 can be described with the following 
function [21] 


Fzero X= 0 
g(x) O<xX<Prate 
(x)= 2 
f Frate X = Prate 2) 
0 otherwise 
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Fig. 2. Distribution of WTG output power. 
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Fig. 3. Multi-dimension distribution of WTGs output power. 
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where function g(x) can either be expressed as a fitted polynomial 
function or be represented by discrete samples. Also, in the case of 
existence of two WTG units in a system, WTGs output distribution 
can be depicted as Fig. 3. 


2.1. Power system analysis in the presence of uncertainties 


In electric power systems installed with WTG, as mentioned in the 
previous section it is necessary to consider WTG uncertainties. Also, 
other kinds of uncertainties such as load uncertainties should be 
incorporated. Here, a general method to tackle with uncertainties is 
presented. The scenario based MCS is a commonly used and suitable 
tool for incorporating uncertainties in problem. The proposed algo- 
rithm as shown in Fig. 4 is a general method which can be used to 
tackle with any kind of uncertainties in power systems. At first, for 
each scenario of MCS, the wind speed and load are randomly 
generated based on their probability distribution. In this paper, the 
uncertainties of the load are modeled as a normal Probability 
Distribution Function (PDF). Then, the WTG output and its probability 
can be calculated by using (1) and (2) respectively. Also, the probability 
of the random load can be achieved by using its normal PDF. In this 
situation, the WTG and load are set on their randomly generated 
values. 

By using the load probability and WTG output probability, the 
probability of the current scenario is also calculated. Then, an 
Optimal Power Flow (OPF) is carried out and the required outputs 


y 
Generating wind speed and load 
value based on their distribution 


á d 
Calculating the WTG output based on the WTG 
input-output relationship depicted in Fig.1 
v 


Calculating the WTG output probability by using 
equation (2) and load probability by using its PDF 


| 


Calculating the probability of scenario by using 
WTG output probability and load probability 


Vv 
Performing OPF 


y 
Recording the required varaiables such as line 
transmission powers, nodal prices, social welfare 


Converge? 


Calculating the PDF function for each variable 


Fig. 4. Flowchart of generating PDF by using MCS. 


such as line transmission powers, nodal prices, and social welfare 
are recorded. Thus, in each scenario, the variables and their 
probabilities are obtained and PDF can be easily calculated. The 
scenarios are iterated until a breaking criterion is met. By using the 
proposed method, the relationship between WTG output powers, 
the probability of WTG output powers and the system parameters 
are obtained. Fig. 5 shows the relationship between WTG output 
power and system congestion surplus in test case 1 (test systems 
are presented in section 5). Also, Fig. 6 shows the relationship 
between WTG output power and line 2-3 transmission power in 
test case 1. It is clearly seen that the WTG uncertainties affect the 
system parameters and there is a PDF for each parameter of 
system such as line transmission power or system congestion 
surplus. With regard to this issue, Power system planning in the 
presence of PDF for parameters is very complicated. 

Fig. 7 shows the relationship between WTGs output power and 
system congestion surplus in test case 2. In this test case, there are 
two WTG units. It is seen from Figs. 5-7 that the WTG output has a 
great effect on the system parameters and this uncertainty should 
be incorporated in any planning. Thus, it is valuable to find a way 


500 + 4 
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300 + 4 
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Congestion surplus ($/MWh) 


100 + 4 


-20 0 20 40 60 80 100 120 140 160 
WTG output power (MW) 


Fig. 5. Relationship between WTG output power and congestion surplus in test 
case 1. 


Transmission power form bus 2 to bus 3 (MW) 


-20 0 20 40 60 80 100 120 140 160 
WTG output power (MW) 


Fig. 6. Relationship between WTG output power and transmission power of line 
2-3 in test case 1. 
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Fig. 7. Relationship between WTGs output power and congestion surplus in test 
case 2. 


to convert the PDF function of a parameter to an equivalent fixed 
number. The next section presents a mathematical-probabilistic 
method to deal with this problem. It is worth mentioning that, 
Figs. 5-7 show the relationship between WTG output power and 
system parameters. But, Figs. 2 and 3 show the relationship 
between WTG output power and the probability. So, the relation- 
ship between system parameters and their probability (this rela- 
tion is known as PDF) can be easily calculated. In fact, there is a 
correlation between WTG output power, the probability of WTG 
output power and any system parameters. 


2.2. Calculation of expected value for parameters with PDF 


Expected value (mean) is a very effective and commonly used 
method to convert a PDF to a number. It is the weighted average of all 
possible values of a random probabilistic variable. By calculating the 
expected value, not only the nature of uncertainty is considered, but 
also the calculations and simulations can be performed without extra 
complexity. The expected value of the variable X which takes value x, 
with probability p; is calculated as follows: 


n 
E(X)= X XkPk (3) 
k=1 
Therefore, the system uncertainties can be eliminated and any PDF is 
converted to a number. The proposed method is used to include the 
parameters with PDF (such as market congestion surplus) into TEP 
problem with considering the wind power uncertainties. 


3. Problem formulation 


In electric power systems, technical and economical aspects are in 
conflict with each other. Therefore, considering one of them as 
objective function is not appropriate. Considering different objectives 
for TEP problem covers all aspects of planning. Thus, in this paper, 
both are considered. Technical objectives try to improve the power 
system operation, but economical objectives try to reduce the cost. In 
this paper, the reduction of the investment cost of lines and reactive 
sources are considered as economical objectives. Also, increasing 
social welfare of the network is considered as a technical objective 
function. These objective functions are defined as; 


Cost of TEP problem (vo): the cost on new installed lines is 
computed as (4) [17]. 


vo=c'n (4) 


Cost of RPP problem (vı): the cost of new locally reactive 
sources is computed as (5) [17]. 
vı =f1(q4,u)= $, (Cont Crd Ux (5) 


keQy 


Cost of social welfare (v3): the social welfare of the system is 
obtained by (6) [25]. 


v2 =fo(Ps, Pr, Hi, Us) = XHn x Pi- LTT si. x Psi) (6) 


The objective function və shows the difference between the 
payment of consumers and the revenue of producers. This objec- 
tive function is known as “congestion surplus” [25]. Besides, the 
minimization of congestion surplus is equal to the maximization of 
the social welfare and the two expressions are equivalent in 
meanings. Thus, this paper tries to maximize the social welfare. 

The objectives vo and v4 are per “$”, but v2 is per “$/h”. Thus, 
these units should be converted to the same unit. For easy 
comparison, all parameters are converted to “$/year”. Therefore, 
the investment cost (vo and v1) is converted to “$/year” by (7) [26]. 
Also, the congestion surplus is converted to “$/year” by (8). 


rxa+r" $ 


Vannual = Vtotal X (+n™—1 year (7) 
v zyx 34x365 2 (8) 
2 annual = ¥2 year 


where, the coefficients “24” and “365” show the “hours of a day” 
and the “days of a year”, respectively. Eventually, the objective 
functions are rearranged as follows: 


rx (A1t+nry™ $ 


Jo=vo x (1+n™-1 year (9) 
rex(1+r™ $ 
J=% x (1+n™—1 year a0 
$ 
Ja=¥a x 24x 365 (11) 


where the variable “year” shows 365 days, Jo demonstrates the 
annual cost of new installed lines, Jı demonstrates the annual cost 
of new locally reactive sources and Jz is the annual cost of 
congestion surplus of the system. These three parameters are 
considered as the objective function for TEP in the following form: 


Min Jo (12) 
Min J, (13) 
Min J, (14) 
Subject to, 
PV, O, n)—P¢+Pp=0 (15) 
QV, O,n)—Qc¢+Qp-q=0 (16) 
PE™ < Po < PE* (17) 
Qa" < Qc < QG” (18) 
y™in <V< y™ax (19) 
(N+NO)S!"°" <(N+NO0)S™ (20) 
(N+NO)S" < (N+NO)S™™ (21) 


qn <q< qr™ (22) 
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Ozanan" (23) 


EENS < EENS” (24) 


Constraints (15) and (16) show the power balance in buses. The 
bounds of generator output power are shown by (17) and (18). 
Constraint (19) shows the limits of voltage levels. The limits of 
transmission power in lines are presented by (20) and (21). 
Constraint (22) shows the limit for installation of reactive sources. 
Constraint (23) represents the limit of the new installed lines in 
each corridor and the reliability constraint is represented by (24). 
The parameters P(V,Q,n) and Q(V,@,n) in (15) and (16) are 
computed as (25) and (26). Also, the bus admittance matrices (G 


and B) are calculated as (27) and (28). 

P(V,O,n)=V; 2 Vi[G;(n) cos 0+ By(n) sin Oj] (25) 
jeNg 

Q:(V,O,n)= Vi 2 V,[Gj(n) sin 0+ By(n) cos Oj] (26) 
jeNg 


Gy(n) = —(njZy+ 7929) 
G=9 Gins 2 (nijgj+ nggi) 27 
Jesa 


0 
By(n) = — (nbi + nġbi) 


B=) Bin) abi + 3 naby bi +915 + OON 2) 
jeQ 
Eventually, the element (ij) of S*°" and S® in (20) and (21) are 


calculated as (29) and (30). 


= ee or +Q) 29) 
Sp = y (PP)? +Q) (30) 


where, 

a =Vigi—-ViVi(gq cos 0; +by sin 0) (31) 
ae —V? (bi + by) -ViVi(Bi sin 0;—by cos Oj) (32) 
P? = Vi g,—-ViVi(gy cos Oj +by sin Oj) (33) 
Qi? =-V; (D5" + by)— ViV,(gi sin Oj—by cos 6y) (34) 


4. The proposed solution 


In this paper, a multi-stage algorithm is provided for solving 
the problem and its flowchart is shown in Fig. 8. The solution is 
performed in three stages. In the first stage, the TEP problem is 
performed by using the PSO method. In the second stage, the RPP 
is performed and in the third stage, the reliability of planning is 
checked. 


4.1. First stage: solving the TEP problem 


In this stage, the TEP problem is carried out by using PSO 
algorithm. Blocks A to H in Fig. 8 show the details of this stage. In 
block A, all reactive demands are procured locally. In other words, 
the reactive demands are not considered in this phase of planning. 
Then, in blocks B to H, the TEP problem is performed by using PSO. 
In block B, an initial population is generated randomly. In block C, 
the PDF function of congestion surplus is calculated for each 
particle in the population. The PDF function is calculated based 
on the method presented in Fig. 4. It means that, in block C, the 
flowchart shown in Fig. 4 is performed one time for each particle. 


[ Yy 
A | Assuming all reactive demand sare locally procured 
v 
B Generating the initial population a 
v 
Obtaining the PDF of congestion surplus for each 
C individual by using the flowchart in Fig. 4 
v 
D Calculation of expected value of congestion surplus 
for each individual 
Stage 1 
8 Vv 
E Calculation of (Jọ+J2) as investment cost for each 
individual 
Vv 
F choosing the best individual with minimum cost 
A A 
G Updating the population 
H Converge? ° 
= Y Yes 
lı Installation of the new determined lines 
Y 
I Removing local reactive sources and considering 
Stage 2 2 reactive demands 
Vv 
l minimizing J, by reactive power planning 
J EENS > EENS max? Yes 
Stage 3 LOLE > LOLE max? 
No 


: 


Fig. 8. Flowchart of the proposed planning. 


In a population with n particle, the flowchart in Fig. 4 is performed 
n times and n PDF functions are obtained. Then, in block D the 
expected value of each PDF function of the congestion surplus is 
calculated. The investment cost of each particle is computed in 
block E. Then, in block F, the elite particle with the minimum cost 
(the sum of investment cost of new lines and the social welfare 
cost) is chosen. In block G, the PSO population is updated based on 
the PSO rules and in block H, the stopping criterion is checked 
which is 15 iterations without changing the best elite. The PSO 
algorithm is repeated till convergence. After convergence, the new 
lines are determined. As referred above, at first, all reactive 
demands are procured locally and hence, the cost of RPP ( J4) is 
not considered in this stage. Therefore, the TEP problem in this 
stage is carried out with two objective functions. The objectives 
are the minimization of the investment cost of new lines ( Jo) and 
the minimization of the congestion surplus ( J2). The sum of these 
two objectives is considered as the final objective function in this 
stage. Thus, the proposed TEP problem is formulated as follows: 


Min (Wo Jo+W2 J2) (35) 


Subject to 
Eqs. (15)-(21) and (23) 
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The objective function (35) shows a multi-objective optimiza- 
tion problem, where coefficients wọ and wz represent the weight- 
ing factors. In this paper, Jo and Jz are considered with equal 
importance and wọ=w2. It should be noted that the weighting 
factors are chosen based on the importance of parameters and also 
the system conditions. Since both Jo and J2 have same unit “$/year”, 
their importance is equal and the same weights are considered 
for them. 


4.2. Second stage: solving the RPP problem 


Block I in Fig. 8 represents this stage in details in which RPP is 
performed by using the PSO method. At first, in block l4, the 
obtained new lines from stage 1 are installed. Then, in block In, the 
reactive demands are added to the network. Eventually in block Is, 
RPP is performed by using PSO and the feasible locally reactive 
sources are determined. The RPP is carried out to determine the 
place of locally reactive sources with the minimum cost. In this 
paper, RPP is performed to minimize the investment cost of locally 
reactive sources (J4) and is formulated as the following optimiza- 
tion problem: 


Min J, (36) 


Subject to 

Eqs. (15)-(22) 

At first, in blocks B to H, the TEP problem is solved and then, in 
block I3, the RPP is performed by using PSO. 

It should be noted that the output of the wind power is 
uncertain; therefore the RPP should be incorporated associated 
with the wind power uncertainties. To deal with this problem and 
considering the wind power uncertainties in RPP, the presented 
planning in (36) is performed for the different values of the wind 
output. The RPP is performed for different values of the wind 
output which changes from zero to its nominal value. In the 
proposed method, the output of the wind power changes with 
1 MW discrete steps. Eventually, the minimal reactive sources are 
chosen for installation as locally reactive sources. Fig. 9 shows the 
RPP output versus the output of the wind power in test case 1. It is 
seen that the minimum amount of locally reactive sources are 
obtained, when the output of the wind power is at its maximum 
rate. Therefore, for any RPP, the output of the wind power is fixed 
at the maximum power rate and then the RPP is performed. 


120 T T T T T T T T 
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20 F 


WTG power (MW) 


Fig. 9. Locally reactive sources versus WTG output power in test case 1. 


4.3. Third stage: reliability assessment 


After installation of new lines and new locally reactive sources, 
the reliability of the system is evaluated in this stage. The resulted 
planning with EENS or LOLE greater than EENS™™ and LOLE™®* is 
eliminated and the algorithm is repeated from the beginning. In 
this paper, in order to assess the reliability, a constraint is 
considered for the planning and it is checked in block J. In this 
block, the reliability of the planning is evaluated by using the 
reliability indexes EENS and LOLE, which are calculated using (37) 
and (38). 


D Total amount of loss of load 
Ns 


EENS = (MWh/year) (37) 


E Total hours of loss of load 


LOLE = N; 


(h/Year) (38) 


In order to assess the reliability, the combination of OPF and 
scenario-based MCS is used. In all scenarios of the MCS, the load 
and WTG output are randomly generated based on their estimated 
PDF and also a random number is considered for each line outage. 
For each scenario, the AC-OPF is performed to determine the 
unsupplied power. The EENS is calculated as the average amount 
of unsupplied powers for all scenarios. 

In this paper, PSO and OPF are used in the planning. Although 
these two routines are well known, in the next sections, they are 
briefly explained. The reliability with consideration of DGs is 
reported in [27]. 


4.4. Optimal power flow 


OPF is a well known tool in power system analysis. It is a 
nonlinear programming with equality and inequality constraints. 
It is carried out with different objective functions such as the 
maximization of the distance to collapse, the maximization of the 
social welfare etc. Among these objectives, the social welfare 
maximization is a market oriented objective function. In this 
paper, a market based OPF is used, which aims at the social 
welfare maximization (or the congestion surplus minimization). 
The OPF with the proposed objective function is formulated as 
follows [28]: 


Min Jp (39) 


Subject to 
Eqs. (15)-(21) 


4.5. Updating the population in PSO 


As referred before, PSO method is used twice in the proposed 
planning; at the first stage of the planning, TEP is solved using 
PSO. In this stage, PSO finds new lines to expansion. In the second 
stage, RPP is solved using PSO, where the algorithm finds reactive 
power sources. PSO is a well-known optimization technique which 
has been widely used to solve different constrained optimization 
programming. The algorithm has been inspired from social beha- 
vior of bird flocking. PSO begins with a random population matrix 
and it has no evolution operators such as crossover and mutation. 
The rows in the matrix are called particles. They contain the 
variable values and are not binary encoded. Each particle moves 
toward the cost surface with a velocity [29]. In PSO algorithm, each 
particle is defined as a multi-dimensional particle with two values 
of pig(t) and vjg(t). In each stage of the movement of the swarms, 
each particle with two best values is updated. Also, Ppes:(t) and 
Zbest(t) are the local and global best solutions. After finding Pyes(t) 
and Zpest(t), the particles update their velocity and position as 
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given by (40) and (41); 


Viat +1) = W()Vja(O) + C1 ANA (Presti, a(t)—Pia(t)) + Catan (Spesrq(t)—Xia(t)) 
(40) 


Pia(t + 1) = Pia(t)+ Via(t + 1) (41) 


where rand is a random value in the range [0,1], the inertia factor 
w is used to control the impact of previous velocities on the 
current velocity during optimization process. The factor w is 
linearly decreased from 0.95 to 0.2 [29]. The particles of PSO 
algorithm usually quickly converge at the initial iterations. 


5. Illustrative test cases 


To show the ability of the proposed methodology, two standard 
test cases are considered in this paper; the modified Garver 6-bus 
and the IEEE 24-bus system (RTS). The first test case is a small 
power system, but the second test case (IEEE 24 bus test system) is 
a standard and large scale power system provided by IEEE which 
has been used in many papers as a test system. Generally, the 
proposed methodology does not have any limitations and can be 
used for any practical and large scale power systems. 


5.1. Test case 1 


A suitable system for TEP studies is Garver 6-bus system [2]. 
Sine in the system bus 6 is not connected to the rest of the system, 
it should be connected via new lines. Thus, TEP problem should be 
incorporated. The basic configuration of this system is depicted in 
Fig. 10. This system has five buses, six installed lines and 15 
candidate lines; its total demands are 760 MW and 152 MVAr 
respectively. In market environment, this system has 10 generating 
units and five loads [16]. The total generation power at bus 3 is 
300 MW, and in this paper, half of this power is considered as 
WTG with wind speed parameter values: Va=4, Vrare=10 and 
V«=22 m/s. The wind speed has a mean 5.4 m/s and standard 
deviation 2.7 m/s [21]. It is assumed that the maximum lines 
between two nodes are five. Also, the fixed and variable costs of 
reactive sources are Co = 1000$ and c1 =3 $/kvar, respectively [3]. 
The EENS™ is limited by 1%. The project lifetime and the discount 
rate are considered as 15 years and 10% respectively. Also, the 
maximum and the minimum voltage levels are 105% and 95%, 
respectively. Normal distribution with 10% standard deviation is 
assumed for loads. The FOR and MTTR for all lines are also fixed as 
0.03 (fail/year) and 50 (h/year), respectively. 


Ds 


Fig. 10. Modified Graver 6-bus test system. 


5.2. Test case 2 


In order to show the effectiveness of the proposed planning in 
large scale power systems, the IEEE 24-bus Reliability Test System 
(RTS) is also used. The system is a well known test system with 24 
buses and 41 lines. The system data for power flow and market 
studies are presented in [4], respectively. Two 180 MW WTG are 
assumed in buses 1 and 15. For both WTGs, the wind speed 
parameter values are Va=4, Vrate=10 and V.9.-=22 m/s. Also, wind 
speed has a mean 5.42 m/s and standard deviation 3.06 m/s in bus 
1 and the mean value 5.41 m/s and standard deviation 2.7 m/s in 
bus 15 [21]. The other parameters are considered the same as first 
test case. It is worth mentioning that all electricity market data 
have been obtained from the realistic cases of the day-ahead 
electricity market of mainland Spain based on demand patterns 
corresponding to year 2004 and all data have been matched with 
actual values observed in actual markets [4]. 


6. Results and discussions 


In this section, the proposed planning is carried out on two test 
cases. The proposed constrained optimization planning is solved 
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Fig. 11. Convergence of PSO and GA for the first stage of the planning in test case 1. 


Exiting lines 


=== == New lines 


Fig. 12. New determined lines for test case 1. 
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by using PSO. In order to show the accuracy and the effectiveness 
of the PSO method, the proposed optimization problem is also 
solved by GA and results are compared and studied. 


6.1. Results of test case 1 


The proposed planning is carried out on test case 1. The conver- 
gence of PSO and GA in solving the first stage of the problem is 
depicted in Fig. 11. The convergence is obtained after 15 iterations. The 
system after installation of new lines is depicted in Fig. 12. It can be 
seen that five new lines are installed. The results of planning are 
presented in Table 1. The table shows that, PSO-based planning leads 
to five new lines, while the GA-based planning gives six new lines. In 
RPP stage, 96.3 MVAr is obtained in PSO method, but GA method leads 
to 85.01 MVAr. This means that the GA method leads to less locally 


Table 1 
Results of planning for test case 1. 


PSO results 


reactive sources than the PSO method. Besides, the congestion surplus 
in PSO case is 5% greater than GA case. Thus, PSO method leads to 
fewer new lines, but locally reactive sources and congestion surplus 
are more than GA method. A complete comparison between GA and 
PSO is presented in Table 2. This table shows that, the investment cost 
of new lines in PSO method is 2.6294 ($ x 10°/year) less than GA case. 
But, the investment cost of locally reactive sources in PSO case is 
0.0044 ($x 10°/year) greater than the GA method. Besides, the 
congestion surplus of the PSO method is 0.2268 ($ x 10°/year) greater 
than the GA method. Eventually, the total cost of PSO method is 
2.3983 ($ x 10°/year) less than GA method. Regarding all results, it can 
be concluded that PSO method is a strong optimization technique and 
can solve the proposed constrained optimization problem better than 
GA method. PSO method gives a better planning with minimum cost 
and the obtained plan leads to lower cost while all system constraints 


GA results 


Lines addition in AC-TEP phase 


n2-6=2, N4-6 


2, ns-6=1 Nı-5=1, n2-6=2, N4-6=2, Ns-6=1 


VAR source allocation phase 


Expected value of congestion surplus 


Table 2 
Comparison of PSO and GA for test case 1. 


38.36 MVAr at bus 4 
57.94 MVAr at bus 5 
525.8423 ($/h) 

4.6064 x 10° ($/year) 


41.68 MVAr at bus 4 
43.33 MVAr at bus 5 
499.9572 ($/h) 

4.3796 x 10° ($/year) 


Proposed method with PSO Proposed method with GA 


Number of new added lines 5 6 
Investment cost of the new lines ($ x 10°) 181 201 
Investment cost of the new lines ($ x 10°/year) 23.7968 26.4262 
Investment cost of reactive power planning ($ x 10°) 0.2899 0.2560 
Investment cost of reactive power planning ($ x 10°/year) 0.0389 0.0345 
Congestion surplus ($/h) 525.8423 499.9572 
Congestion surplus ($ x 10°/year) 4.6064 4.3796 
Total cost of planning ($ x 10°/year) 28.4412 30.8395 
Table 3 
Results of planning for test case 2. 
PSO results GA results 
Lines addition in AC-TEP phase: Ng-10=2, N7_-g=1, N4-13=1, N6-17= 1 Ne-10=2, N7_-g=1, N11-13= 1, Ni6-17= 1, Ng-12=1 


206.1 MVAr at bus 3 
49.41 MVAr at bus 4 
139.1 MVAr at bus 9 
573.7 MVAr at bus 12 
131.2 MVar at bus 24 
9138 ($/h) 

80.0489 x 10° ($/year) 


VAR source allocation phase: 


Expected value of congestion surplus 


Table 4 
Comparison of PSO and GA for test case 2. 


Number of new added lines 

Investment cost of the new lines ($ x 10°) 

Investment cost of the new lines ($ x 10°/year) 
Investment cost of reactive power planning ($ x 10°) 
Investment cost of reactive power planning ($ x 10°/year) 
Congestion surplus ($/h) 

Congestion surplus ($ x 10°/year) 

Total cost of planning ($ x 10°/year) 


198.71 MVAr at bus 3 
29.11 MVAr at bus 4 
140.2 MVAr at bus 9 
551.8 MVAr at bus 12 
128.3 MVar at bus 24 
8888 ($/h) 

77.8589 x 10° ($/year) 


Proposed method with PSO Proposed method with GA 


5 6 

150 200 
19.7211 26.2948 
3.2991 2.8461 
0.4337 0.3742 
9138 8888 
80.0489 77.8589 
100.2037 104.5278 
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Table 5 
The reliability indexes for both the test cases. 


Test case 1 Test case 2 

PSO GA PSO GA 
EENS (MWh/Year) 136.94 136.04 11158.3 11165.2 
LOLE (h/Year) 15.60 19.51 81.18 85.41 


are satisfied and also the system performance is guaranteed. Besides, 
PSO technique gives one line less than the GA method and this result 
is very favorable for power system operator due to the problems of 
installing new lines. 


6.2. Results of test case 2 


In this section, the proposed planning is carried out on test case 2. 
Table 3 shows the planning results. Five new lines are denoted for 
network expansion in PSO planning, but in GA planning, six new lines 
are denoted for TEP. In the second stage of planning which is RPP, it is 
seen that the PSO-TEP denotes 1099.51 MVAr locally reactive sources, 
which is 51 MVAr more than the GA method. Besides, the congestion 
surplus in PSO case is 250 ($/h) or 2.19 x 10° ($/year) greater than the 
GA case. Thus, like the first test case, PSO planning gives less new lines 
and more locally reactive sources and congestion surplus. Therefore, a 
complete economical comparison of two methods is valuable and 
presented in Table 4. As seen from this table, the investment cost of 
new lines in PSO planning is 6.5737 ($x 10°/year) less than GA 
planning. But in RPP stage, the investment cost of GA is 0.0595 
($ x 10°/year) less than PSO. Also, the congestion surplus of GA 
method is 2.19 ($ x 10°/year) less than PSO. Hence, the total cost of 
planning can be calculated for complete comparison, where, the total 
cost of the PSO method is 4.3241 ($ x 10°/year) less than the GA 
method. Considering all results and comparisons, it is very clear to see 
that PSO technique can be used to solve such constrained optimization 
programming. PSO method gives an optimal plan with less cost and 
better performance. It can also be concluded that both PSO and GA 
methods can solve the problem, but PSO results are more optimal than 
GA. Eventually with regard to the results, the PSO based plan is chosen 
as final plan. 


6.3. Reliability assessment 


In long term planning, the reliability assessment is an inevitable 
parameter. Thus, the reliability of the resulted planning is evaluated by 
using the most commonly used reliability indexes such as LOLE and 
EENS. The MCS is chosen for simulation as a suitable tool to deal with 
uncertainties. Table 5 shows the reliability indexes for both the test 
cases, respectively. As seen from the table, the reliability of GA-plan 
and PSO-plan is acceptable and both plans are reliable. However, from 
the view of the comparison, the PSO technique leads to a plan with 
less EENS and LOLE, because the obtained plan by PSO comprises less 
components and its reliability is more. 


7. Conclusions 


In this paper, considering load and wind uncertainties in TEP 
problem under electricity market was successfully investigated. The 
presented TEP-RPP problem is a multi-objective optimization problem. 
In the proposed planning, the social welfare maximization, the 
investment cost minimization and the reliability improvement were 
simultaneously obtained. The uncertainties were involved by using 
MCS and a PDF was obtained for each of them. Also, the reliability 
assessment was carried out by using MCS. Two optimization methods 


were compared and investigated. The simulation results demonstrate 
that, the PSO method show a significant advantage than the GA 
method. Also, the proposed method significantly improves the social 
welfare of the system and resulted in fewer new lines. 
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